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Abstract 

For a volume preserving map, the exit time, averaged over the incoming set of a region, is 
given by the ratio of the measure of the accessible subset of the region to that of the incoming set. 
This result is primarily of interest to show two things: first it gives a simple bound on the algebraic 
decay exponent of the survival probability. Second, it gives a tool which permits the computation of 
the measure of the accessible set. We use this to find the measure of the bounded orbits for 
Henon's quadratic map. 

1 Introduction 

Orbits of an area preserving map are eternally trapped in a region if they are enclosed by an 
invariant circle, but can leak through destoyed circles (cantori). However, even untrapped orbits feel 
the influence of these invariant boundaries because they are extremely "sticky."[l-13] In higher 
dimensional maps, the mechanisms for trapping and escape are much less understood [14]. In this 
paper we study the time of escape for orbits that begin in a specified region, A. We discuss the exit 
time, the time for a point to first exit the set, and the transit time, the time for a point to traverse the 
set. The exit time distribution, is the probability distribution of exit times. Our primary goal is to 
use this distribution to probe the trapped invariant set. 

The study of escape time for chaotic maps has long been of interest. As an example, 
consider Henon's area preserving, quadratic map, which we write as 

H: (x,y)^(y-k +x\-x) . (1) 

We are interested in the set of bounded orbits. It is possible to show that (for k>-l) all bounded 
orbits are contained in the square { z = (x,y): \x,y\ < Xs} where x,^ = 1 ± \/l + A: • Here we denote the 
two fixed points (for k>-l) by Ze = i^er^e) and Zs = (xs, -Xs)- The first is elliptic when -1<^<3, and 
the second is a saddle. Bounded orbits include all periodic orbits and all orbits within invariant 
circles; the latter exist, for example, in the neighborhood of Ze providing it is elliptic and has a 
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rotation number, CO, which is not V3 or V2 (i.e., k is not ^/4 or larger than 3). There are also trapped 
hj^erbolic orbits, for example those that are homocUnic to the saddle fixed point. 

In his original paper [15], Henon studied the set of bounded orbits of H (in a different 
coordinate system than ours) by iterating points along a segment of the symmetry line {x=-y) to 
determine the subintervals that did not reach a fixed large distance from the origin within 100 
iterates. He noted that as k varies, the boundaries of the trapped intervals either closely follow the 
position of an unstable rotational periodic point or else an island around an elliptic periodic point. 
Channon and Leibowitz [16] studied the escape and trapping from the period 5 island chain in the 
Henon map at A: = -0.422. They identified the exit and entrance lobes of the resonance as the 
important sets to consider. Studying the power of the map, they started 7750 orbits in the outer 
entrance lobe and computed the survival probability distribution. This was found to decay as t'^ 
with a = 0.5 for short times (up to 10 iterates), but deviated from this for moderate times (up to 45 
iterates). 

Kamey [17] also considered the Henon map, and studied the survival probability for a 
square A enclosing all bounded orbits. He mostly studied k= -0.6 where the most prominent island 
chain is period 6. Kamey' s primary object of study was the "trapping time statistic" which is 
proportional to the exit time distribution for A. He found that this distribution decayed as ?~0f+2 
with a about 0.25, for times up to 10^ iterates, though the slope appears to vary and not settle down 
to a fixed value for times up to 10^. Chirikov and Shepelyansky [18] did similar computations for 
the standard map for the parameter value when there is a critical golden circle. They obtained a = 
0.34 up to 10^ iterates. The exponent for the decay is apparently not universal [19]. One reason this 
might be so, according to Murray [20], is that the self-similar limit is not reached for "short" time 
computations. 

Rom-Kedar and Wiggins have emphasized the fact that one can obtain a complete 
description of transport through a region by considering the future history of only the entering 
trajectories [5]. Using their notion of "lobe dynamics" for a homoclinic tangle, they obtained an 
expression for the accessible region in terms of the exit time distribution — we will use a slightly 
generalized version of this below. 

2 Definitions 

Let/;M ^ M be a homeomorphism with an invariant measure |l. For measurable sets A ,B 
cz M, the crossing time (or first passage time) for a eAioB is defined as 




(2) 
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We let tA^Bici-) = if fl never reaches B. If B is the complement of A, then the crossing time is 
called the (forward) exit time 



t+(a) = tA^M\A(a) . 



(3) 



Similarly the backward exit time for a e A is defined as 



t (a) = min 



n>0 



(n:f-"(a)^A) 



(4) 



lfB=A, then the crossing time is called the first return time: 



tremmia) =?a^a(«) . fov a G A. 



(5) 



According to our definition t,^,^r„(a) > 1, and is equal to 1 whenever the point does not leave on the 
first iterate. We define the exit set £ cr A as the subset with exit time 1, and the entry set / c: A as 
the set with backward exit time 1; equivalently. 



The turnstile is the union of E and /. 

As an example, consider the Henon map (1). A natural choice for the region A is the 
"resonance zone", shown in Fig. 1, consisting of the region bounded by the segments of the left- 
going branches of the stable and unstable manifolds of the saddle fixed point up to their first 
intersection on the symmetry line {x=-y). This region contains all bounded orbits of H. 

The transit time of a point a is the sum of its forward and backward exit times minus one (to 
get rid of an annoying term). 



The transit time is an orbit invariant: each point along an orbit has the same transit time. The transit 
time decomposition of a set is the partition of a set into subsets with equal transit time. We denote 
the transit time decomposition of the entry set / by sets 7); this is also the part of I with exit time j: 



E = AYI(A),I = AY(A). 



(6) 



ttransit(a) = t'^(a)+t'(a)-'^ • 



(V) 



Tj={aGl: t\a)^j}. 



(8) 
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Figure 1. Exit and entry sets for the fixed point resonance of the Henon map for k = 0.5. 



For example, we show the first few preimages of the exit set, in sequentially lighter shades, for the 
Henon resonance zone in Fig. 1. For this parameter value Ti = 72 = = (j) and is shown. A 
partial transit time decomposition of the Henon entry set is shown in Fig. 2. 

(color figure attached) 

Figure 2. Exit time decomposition of the entry set of the Henon Trellis at ^=0.5. Color scale is 
given at the bottom. 

By measure preservation, almost every point that enters A must eventually escape, so 



The accessible set, Aacc A is defined to be the set with finite backward exit time; it is the 
set that can be reached from the outside: 



Of course, the set with finite exit time differs from Aacc at most by a set of measure zero, since the 
set that enters but never exits must have measure zero. 



7 = 1 



(9) 




(10) 



3 Average Exit Time 

In this note we are interested in these transport times averaged over sets of initial conditions. 
For a function g(z) we denote 

Remarkably, there are some simple formulae for average transport times. The following lemma was 
stated for volume preserving flows (without its elementary proof ) in [21] and is implicitly obtained 
for two dimensional maps in [22]. 

Average Exit Time Lemma: The average exit (transit) time for incoming orbits is 

(t^=(t .\=i^^^^ (11) 

V // Y transit/ 1 ^(-J-) K^'^J 

Proof: Since 7} cr/ are disjoint, and cover almost all of 7, and the exit (and transit) time of the set Tj 
is j, the average exit time is given by 

(•').= w>% ' ''^'''^ ■ 

(assuming the sum exists). Given the transit time decomposition Tj, we can compute the accessible 
subset of A using the sets 

Tj =fXTj) , 

Note that Tj a A for / = 0,...,j-l, and that these sets are disjoint because they have transit time j and 
backward exit time z+1. Furthermore the union of these sets is, up to measure zero, the entire subset 
of A that exits. Thus by definition 

A,,,=(juTj ^ ii(A„j=|: 1:11(7;) . 

7 = 1 (=0 j = l 1 = ■' 

Since measure is preserved [i{Tj) = ^(Tj) , and so 

[liA^j-tim (12) 

This expression was first obtained by Rom Kedar and Wiggins (see e.g., (6.13) in [22]). Finally 
since fM^acc) ^iM^)-, it is clear that the sum converges. □ 

Since iM^acc) is finite, a simple consequence of this result is 
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Corollary 1: The measure of the region with transit time t must decay at least as t"^. 

Furthermore, the lemma implies a well known result of Kac in ergodic theory [23]. In our context 
this can be generalized to 

Corollary 2 (Kac): Suppose n(M)=l. The average first return time to a region A e M 
is 

t ) (13) 

''return/ A ^(ji^^ ' V^-'-' 

where Mace is the subset of M that is accessible to orbits beginning in A. 

Proof: For points that stay in A for at least one step, the first return time is one. The remaining 
portion of A is its exit set E. No'w, f(E) is the incoming set of M\A, and by our Lemma, the average 
transit time for these points through M\A is 

^ _ |l(M„,,\A) _ H(M„ J - ii(A) 



Since there has been one step already, the first return time to A for these points is one larger. So the 
average first return time to A is 



t ) =^ 

•'return/A jx(A) 



|i(A) - |i(£)) X 1 + |i(£) X 



/(£) 



This reduces to the promised result. □ 

We can also easily compute the transit time averaged over A — since each image of Tj has 
the same transit time, the area of A that has transit time j is jfJ(Tj), therefore 

Corollary 3: The average transit time for points that do escape from A is 



providing that this sum converges. The average exit time for A is 



transit) A ^ 



(15) 
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In contrast to the exit time averaged over /, we cannot use disjointness to show that (14) or (15) 
converges. In fact as we show by a simple example in the next section, these sums need not 
converge. 



4 Exit Time Distributions 

As was emphasized by Rom-Kedar and Wiggins [5], if one knows the Tj, one has most of 
the information one could want about transport though A. Here we recall some definitions of 
normalized exit and transit time distributions. 

The exit time probability distribution for the entry set is the probability that a trajectory in / 
will have exit time t: 

Prob{t\I) = j) = ^ . (16) 

Similarly the survival probability, is the probability that the exit time will be at least than k: 

Prob{ t\I) >k)=-^t |i(r,) (17) 

Note that P{t+{I) > 1) = 1 by eq. (9). 

Once we know the 7y, distributions for A are also known. For example, since /(r^+/) (^Aacc 
has exit time k (with backward exit time j+l), the subset of Aacc that has exit time k is given by 
O f\Tk+ i) • Thus the exit time probability distribution for A^^^c is 

ProH nA^,) = k)^ tttJ-^S m = -^Prob{t\I) > k) 

which is the same as the survival distribution for / (17), up to normalization. Similarly the subset of 

J J, 

A with transit time 7 is LJ / (Tj) , thus the transit time probability is 



Probi t,,„„JA^J = = Z |l(r,.) = j ' 

Finally the survival time distribution for A is 

Note that these equations imply that if , e.g., 

11(7^:) ~ H2+«), ask^oo, 

where coO then we have 
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Prob(t\I) = k)~k-^''^ ^ 

ProbitV) >k)~ Prob{f{A,J = k)~ Prob{t,,,,,(A,J =k) ~ k 
Prob{t\A^,J>k)~k-'' 



5 Examples 

Consider the linear, area preserving, hyperbolic map 



X \ _ 

y 



I 

r' 



where ?i > 1. Let A be the unit square {ix,y) : < x,y < 1} . Then the entrance set is the rectangle / = 
AY(A ) = { (x,y) : 0<x< 1, A'^ <y ^ 1 } • It is easy to see that the transit time decomposition of / is 



Tj = {(x,y) : Pi-J <x< X^-J ,X-l<y<l} 
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Figure 3. Transit time decomposition of the entry region for the linear hyperbolic map. The 
transiting regions Tj, j=l,2,3 and a few of their iterates are shown. 

So the measures of each of these regions are 



a-Dd-A') 



(18) 
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These decay exponentially, as one would expect. The calculations needed for Eqs. (9), (12) and (14) 
are derivatives of simple geometric sums, yielding 

tlLi(Tj) = (l-X-') = [i(I) , 
tj^(T^) = l=[l(A) , 



as required by the lemma. Thus the average transport times are 

^ \ _ X + 1 /^+\ _ X 



These results are unchanged if we scale the size of A (since the map is linear), or if we replace A by 
a square centered on the fixed point (since the map is symmetric). 

More generally, suppose the map is a diagonal hyperbolic matrix with eigenvalues (A,i, A,2..., 
A,d, Hd+b 1^2—, |An)> where the Ai>l and the < 1. Then the entry set of the unit hypercube A is the 
union of rectangular cylinders: 



1 = ^{(xi,...x„) sA: ^i, < X, < l| 



d 



If we define A = FI X, , and n= Yl ^, to be the total expansion and contraction, respectively, then 

i = 1 i = d+l 

a simple calculation gives the measure of the transit regions by 

(A- 1X1 -ID 

in complete accord with Eq. (18). In the volume preserving case, we see that AD = 1, and so the 
formulae Eq. (9), (12) and (14) again hold, and the average times are given by eq. (10) with X 
replaced by A. 

Though the region A that we considered above is special, a theorem of Easton [12] implies 
that the rate of escape for any isolating neighborhood of an invariant set are asymptotically similar. 
Thus any region surrounding the fixed point will have escape that is exponential with rate A. 
Similarly, any diagonalizable hyperbolic map with expansion A will also have the same asymptotic 
rate. 

We expect exponential decay of the transit time decomposition for a hyperbolic system. 
Numerical observations of transport, however, indicate that the transit time decomposition decays 
algebraically when there are elliptic regions in A. The simplest example of this behavior is the trivial 
shear: 
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x =x-\-y 

y' = y 



Let A be the unit square as before. Now the entry set is the triangle / = { {x,y) = 0<x<y< 1 } . The 
transit time decomposition is 

Tj = l^(x,y)e I:l-jy <x <l-(j-l)y| . 



r 




Figure 4. Transit time decomposition for the simple shear. 



These sets have measures which decrease algebraically 



j(j -1) 

The sums to get the average transport times are elementary telescoping sums, and again these verify 
Eqs. (9) and (12). However, for this example the average transit time for A, (14), does not exist. 

6 Bounded Orbits for the Henon Map 

As a final example we use the average transit time, which is straightfoward to compute, as an 
effective method to obtain the accessible area, which is not otherwise computable. Here we do this 
for the resonance zone of the Henon map. The calculation involves several steps. First we find the 
points on the minimizing and minimax homoclinic orbits, Zm and zh, which bound the lobe, see Fig. 
1 . We then construct the boundary of the incoming set, by discretizing and W'\ as graphs y'^(x) 
and y\x), to a resolution h = {xm-Xh)IN for a fixed number of pixels N. Generally, in our calculation 
we used N=\0'^. By reversibility, the exit set is the reflection of the entrance set about x+y=0. The 
average exit time, <?+>/, is given by an integral of the piecewise constant exit time over the incoming 
set. We do this double integral in the most naive way by first integrating T{x) = f{x,y)dy for 
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a fixed x, and then integrating over x using Simpson's rule. To compute T{x), we use bisection to 
zoom in on the discontinuities of f*": first evaluate f*" on a grid of size h; ift^ does not change 
between two grid points, we assume (possibly incorrectly) that it is constant between. If there is a 
change, we bisect the interval until either f*" is equal on the endpoints, or the perceived error in 
neglecting the variation in r*" is small enough (we chose an error of 10"^ for this). Then r(x) is the 
sum of the f*" values times the interval lengths. The resulting average exit time is shown in Fig. 5. 
While <?+>/is generally decreasing, there are numerous small upward jumps. We discuss these 
more below. 




Figure 5. Average exit time for the Henon map as a function of parameter. Most points used 
N=IQ^ and 

^max^ 10^, though for k<-.95, we used ?^a;(;=5xl0^ 

The area of the lobe is either given by summing the number of pixels in the exit set, or taking the 
difference in action between the two homocUnic points [24]. Then \l{Aacc) is obtained from Eq.(l 1). 
It is interesting to compare this with the total area of the resonance zone itself, |x(A). This is most 
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easily computed by taking the difference in action between the action of the minimax homoclinic 
point and the fbced point. These areas are shown in Fig. 6. 




Figure 6. Resonance and Lobe area for the Henon map. Also shown is |a(Ti), which is nonzero 
beyond the formation of the geometric horseshoe at k « 5.706. 

Note that the resonance and lobe area grow monotonically and smoothly. The value k = -I 
corresponds to the saddle node bifurcation, where the fixed point resonance zone is created. 
Slightly above this point, the lobe area, AW is apparently exponentially small and most of the 
resonance is filled with invariant curves. 

Combining these results, using Eq. (11) gives the accessible area. In Fig. 7 we show the 
accessible fraction, |i(Aflcc)/M^(^)- The area that is inaccessible, which is identical to the measure of 
the bounded orbits is given from (1 1) by 

ll(Ai) = IJ(A) - ^Aacc) = I^A) - ^\)<t+>i , 
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This area is shown in Fig 8, which is nearly identical to Fig 4. of [1], but our method allows us to 
compute the results to much higher accuracy. The accuracy can be seen better in the next figure, 
which is the best representation of this information. Fig. 9. This shows the inaccessible fraction. 




-10 12 3 4 

k 



Figure 7. Accessible fraction for the fixed point resonance of the Henon map. 

In Fig. 9, the cut-off at a relative area of 10"^ is an artifact of our numerical method — it is quite 
difficult to reduce the error significantly. To do this one must increase the maximum number of 
iterations, tmax, to pick out narrow chaotic layers near invariant circles, and one must increase A'^ to 
find all possible discontinuities in the exit time. It would be nice improve the accuracy by detecting 
these discontinuities and find a form that fits T(x) in their neighborhood, but we have not done that. 
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- 1 1 2 3 4 

k 

Fig. 8. Measure of the bounded orbits for the Henon map. 

The most prominent features in Fig. 9 are local minima near k = and k = 1.25. These 
correspond to the quadrupling and tripling bifurcations of the elliptic fixed point. In the period 4 
case, a pair of orbits with rotation number V4 are created at the bifurcation. Interestingly for 
0<k<OA, the saddle period 4 orbit, which has points on arranged on the square with comers 
±Vfc), has nearly coincident stable and unstable manifolds — for practical purposes, a saddle 
connection. Fig. 10. Furthermore for 0.2<k<0A this feature dominates the inaccessible set. Using 
this approximation gives A/ ~ 4k for A:>0 — this approximation is shown in Fig 9, and gives good 
agreement with the computed results. Of course, there are other inaccessible islands, most 
importantly islands around the elliptic 03 = V4 orbit. These cause the fraction for 0<^<0.2 to deviate 
from our simple form. 
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Fig 9. Inaccessible fraction of the resonance zone for the Henon map. The two solid curves 
represent simple approximations to the area enclosed by the period 3 and period 4 saddle orbits. 
The fraction falls off on average exponentially with k, aside from dips near prominent bifurcation 
points. Along the top are shown bifurcations of the elliptic point (labeled by rotation number p/q), 
and homoclinic bifurcations corresponding to the creation of type 2, 3, and 4 trellises. 

A pair of period three orbits is created by saddle node at k=l. At k=^/4. the period three 
saddle collides with the elliptic fixed point and there are no encircling invariant curves. Near this 
bifurcation, the most important feature is the virtually perfect saddle connection of the manifolds of 
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the period three saddle (see Fig. 10), which has points at ( -(3, |3) — > ( -1 + P, P) — > (-[3,1 - P) , with 
P = V'^- 1 . Using a triangle as an approximation for this area gives Ai ~ V2 (2P-1)2. The resulting 
curve is also shown in Fig. 9; it fits remarkably well for 1.28< k< 2. Note that, contrary to the 
impression given by Fig. 9, the inaccessible fraction does not go to zero at A^l.25. In particular, 
there is a second period three orbit, (V2,-V2) {-^li-^li) (V2, ^li), that is elliptic. The island 
chain surrounding this orbit, see Fig. 1 1^ has a relative area about 5(10)"^, accounting for most of 
the inaccessible area. There are no other islands visible. 




Fig. 10. Henon map for k = 0.2 and k= 1.6. The outermost invariant curve surrounding the elliptic 
fixed point is closely approximated by a saddle connection of the co = V4 and the V3 saddle 
orbits, respectively. Bounds for the two figures are (-1,1) x (-1,1) and (-1,0) x (0,1), respectively 



There are a number of other sharp drops in the inaccessible fraction, the most prominent 
occur at A: = -0.61, -0.413, 0.585 , 0.78, 2.50, and 3.17. These occur when an invariant circle is 
destroyed, suddenly opening up a new accessible region. The newly opened region will be "large" 
if the critical circle is just outside a large island chain. For example for k = -0.414 there is an 
invariant circle just outside the CO = V5 island chain. This invariant circle is destroyed by A: = -0.413, 
leading to the opening of a new accessible domain, and consequent decrease in measure of the 
bounded orbits. 

7 Conclusions 

We have shown that the average exit time from a region is given exactly by the ratio of the 
area of the accessible portion of the region to the area of the exit set. (1 1). It is interesting that this 
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provides a justification for the oft used estimate that an "escape rate" from a region is given by the 
inverse of this ratio. 




Figure 11 One of the islands around a period 3 orbit at A: = 1.25, the period tripling point of the 
Henon map. Plot bounds are (-1.51,-.597) x (-1.48,-.375). 

We use this to provide a nice numerical tool for computing the measure of the bounded orbits for 
the Henon map. Unfortunately, computational resources limit our accuracy in this calculation to a 
relative measure of about 10"-^. 

For the future, it would be nice to apply (11) to study the bounded orbits for higher 
dimensional maps, for example Moser's canonical form for the quadratic symplectic map [25]. 
Such maps are important in appUcation to particle accelerators. 
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